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TECHNICAL NOTE 2451 


MATHEMATICAL IMPROVEMENT OF METHOD FOR COMPUTING POISSON 
INTEGRALS INVOLVED IN DETERMINATION OF VELOCITY 
DISTRIBUTION ON AIRFOILS 
By I. Flugge-Lotz 


SUMMARY 


The Poisson integral involved in the determination of the change 
in velocity distribution resulting from a change in airfoil profile in 
parallel incompressible flow is solved. 

First, three well-developed numerical methods of evaluating this 
integral, all based on the division of the range of integration into 
small equal intervals, and the difficulties involved in each method, 
are discussed. Then a new method, based on the use of unequal intervals, 
is developed, and compared with the other methods by means of several 
examples. The new method is found to give good results for both the 
direct and inverse airfoil problems and is easily adaptable to rather 
complicated problems. It is particularly recommended for all those 
functions where steep slopes in small portions of the region, to be 
integrated exist. 


INTRODUCTION 


The ordinary thin airfoil at small angles of attack produces only 
slight disturbances in the flow of a parallel incompressible fluid. 
Hence, the influences of camber and thickness upon the velocity distri- 
bution may be computed independently and their effects superimposed. 

The effect of camber may be represented by vortices distributed along 
the chord line of the airfoil section; the effect of the thickness, by 
sources and sinks also along the same chord line. The velocities pro- 
duced by these singularity distributions enable one to compute the 
pressure distribution on the airfoil rather quickly. 

Allen (reference l) has presented this singularity method in a 
form which has proved to be very practical for common usage. However, 
in special cases the unavoidable evaluation of the Poisson integral in 
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the course of the computations has given rise to numerical difficulties. 
Such integrals are usually computed "by the application of finite differ- 
ences using intervals of equal length . However, changes in airfoil 
shape, which result in marked changes in the function to he integrated 
in only sms 1 1 portions of the range of integration, require that extremely 
small interval sizes he employed in this range, and, consequently over 
the entire range of integration. This leads to a considerable amount 
of computatio nal work; hence, it appears reasonable to discuss the 
possibility of employing intervals of varying lengths for the evaluation 
of the Poisson integral. 

This investigation was prompted by the difficulties arising from 
the problem of small c hang es in the shape of symmetrical airfoils at 
the angle of zero lift. The examples included in the present report 
are restricted to this case, but the results obtained are in no way 
specialized and may be applied to all problems wherein the PoisBon 
integral occurs. 

This work was done at Stanford University under the sponsorship 
and with the financial assistance of the National Advisory Committee 
for Aeronautics. 

The author wishes to express her appreciation to Mr. H. Norman 
Abramson for his intelligent and skillful help in the computational 
work and for his assistance in writing the final report. The author 
also wishes to extend her thanks to Mr. R. E. Dannenberg and the 
computing staff of the 7- by 10- foot wind-tunnel section of the Ames 
Aeronautical Laboratory, Moffett Field, California, for preparing the 
extended tables of the functions j no and j no *. 


DISCUSSION OF PROBLEM 


The basic reference profile may be given by y^- r = f(x), and its 

velocity distribution may be known from an earlier computation. The 
problem at hand is that of determining the change in the velocity dis- 
tribution resulting from a change in the shape of the profile (indicated 
by the dotted line in the following fig.). The difference of these two 
shapes is designated as Ay^.. 
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Allen (reference 1, p. 7) gives for the change of velocity the equation 


M*o) 

Vo 



( 1 ) 


where V 0 is the velocity of the basic parallel flow. If, by 
conformal mapping of the outside flow region, the center line of the 
profile is transformed into a circle by the relation 


x = | (1 - cos 0) (2) 

then the profile is transformed into a curve approximating the circle 
shown below. 
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The change in velocity due to a change in form will then he given as 


defining 


AT fo\ 

1 


^ cot 

0 - 

Vo ( e °) 

“ " 2 * 

u 

1 0 

dx 

2 


4 ( Ay t) 





dx 


dx 



— — 

it +0 

— — 

it-© 


d0 


( 3 ) 


This is the form most often used for computation purposes^because the 
inverse problem (that of computing the change in shape due to a change 
in velocity dis tribution) utilizes the analytic form 


d (Ar-fc) 

dx 


*o 





did 




defining 


• (r) - (r) 

VO'rt+0 VO/jr-0 

which is strikingly similar. The corresponding formula in the original 
x,y coordinate system is given by 


dx 



Av 

Vp 

X - x 0 


\fa( c - *0) 

)/k(c - x) 


dx 


( 5 ) 


The evaluation of equation (3) may be accomplished by any one of 
several different methods; however, all of these methods employ the 
device of replacing the integral over the range 0 to 2 it by a sum of 
integrals over intervals of equal length AS. The equally distributed 
po in ts 9 n have corresponding values x n which are not equally 

distributed (see following fig.). 
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This arrangement is sometimes favorable 

d ( A yt) 

dx 


upon the particular form of 


, and sometimes not, depending 
(See discussion following 


equation (43).) 


functions 


The use of the angular coordinate 0 has the advantage that the 

or Az are periodic functions in 9, and this 
dx V 0 


periodicity facilitates the organization of the numerical computations. 
The disadvantage arises from the fact that these functions are usually 
given as functions, of x, and, since the analytic form is not usually 
known, any transformations made will lead to small errors. For example, 

if d (fr^ t ) . or AX. f 8 given at special points which do not correspond 
dx ' ' ' 

to 0 n = nA0, then the computor must obtain the values of these functions 
for the values 0f, 0g, an( i so forth by interpolation. 


DISCISSION OF SOME OF THE EXISTING NUMERICAL 
SOLUTIONS OF POISSON INTEGRAL 


The difficulty encountered in the solution of the Poisson integral 

( 9 - JA - 1 


arises from the fact that the term cot 


or 


x - x r 


(equa- 


tions (l) and ( 3 ) , e.g. ) approaches infinity wnen 0 approaches 0 or 


when 
the function 


approaches 

^Art) 


dx 


or 


The difficulty is of much less consequence when 
A v 


is given analytically than when a numerical 


computation is undertaken. As a consequence, any simple integration, 
performed by replacing the integral with a summation over smaller intervals, 
always requires that the interval in which 0 O or x Q is located be given 

special consideration (see following fig.). 
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H 1 H 


^m 




* — I 1 H- 


3: - 




*nH-l 


H 1- 


0*/ 


V, 


4 1 h 


s m+l 


A majority of the solutions currently in use have "been developed 

to such an extent that, for example, ~(^o) I s given by a sum of 

f d(Ay+V\ * ° 

products of single values of 1 ~ ~j and known factors A^j that is, 

/ n 


r ( 0 °) ■ 

* r\ 


i r d fei cot a9 


2* 


10 


dx 


1 2it-0 


2ir 


° d ( A yt) 9* * 

— — fi- cot — d0 
dx 2 


2rt Z_, 




0+1 5^1 cot £ afl * 

t- dx 2 


n 


(6) 


which leads to (see reference 2, e.g. ) 




K^t) ' 


dx 


( 7 ) 
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The coefficients A no depend upon the particular method of 

a fe! 


numerical integration which is employed. If, for example. 


dx 


is 


replaced by a step-curve, that is, assumed constant in every interval 
(see fig. below), one set of values of A RO would be obtained. 



Greater accuracy wbuld be obtained by the assumption that 


d (Ay t ) 

dx 


is 


replaced by straight-line segments (see fig. below), in which case a 
second set of values of A QO would be obtained. 



d (Ay + ) 

A further refinement would be that of assuming — » — to be composed 

dx 

of segments of parabolas, and so forth. Since the accuracy of the 

resulting values of ~ depends upon both the character of the approxi- 

’o 

mate curve and the size of interval taken, it is apparent that the same 
degree of accuracy might be achieved from many different combinations 

of interval sizes and approximations to the function 


d (Ar t ) 

dx 
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Ka lman (reference 2) has used Simpson’s rule for computing the 
Poisson integral , which corresponds to the replacement of the product 


cot 


6 


( 


d ( A yt) 

L dx 2_ 

i.e., where cot 


* 

5 


by segments of parabolas. The "critical interval" 


n* \ 

— col vas carefully treated by using differences 

2 d ( Ay t) N 

of hitler order I incl udin g the fifth derivative of 


V 


dx 


Kalman 


divided the period of 2it in 20 or 40 intervals and calculated the 
corresponding sets of values of A aQ . Other workers at the KACA have 

extended the calculation of these values of A QO to 80 and l60 intervals 
(unpublished information) . d 

Obviously, the time required for computing increases with the 

Vo 

number of intervals taken because of the increased number of multipli- 
cations to be performed. In addition, greater preparations for the 
computing proce oa are necessarily involyed, particularly since the 


values of 


dx 


needed must usually be obtained by interpolation. 


This interpolation has to be done rather carefully as it is often not 
sufficient simply to take the values of the plotted curve of 

d ( Ay t) 

. This curve should be checked by- difference tables if the 


dx 

values 




dx 


are to represent a smooth curve. 


Jn 


For those functions of Ayt which may be well-represented by a 
Fourier series, there exists a sinqple method of evaluating the Poisson 
integral which has apparently been overlooked until the present time. 
This method has the advantage of leading to a computation which does 
not involve the derivative of Ay^-. 


d Kaiman has also suggested a second method for computing the 
Poisson integral (see reference 3). In this second method he uses 

d (Ay t ) 

dx 


Fourier polynomials to represent the function 


The computing 


procedure is very simple; however, the results depend largely on the 


degree of approximation of 


dx 


by such a polynomial. Thus, for 


large families of functions results are good; however, cases are known 
to the author where results were not satisfactory because regions with 
steep gradients may not be represented well enough by a Fourier 
polyn omial of moderate order. 
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Equation (6) may be written in a different form (reference 1, 
equation (43)) as follows: 


Av = 1 r d ( Ay t) sin _g 
V 0 jrjo dx cos 0 - cos 0 Q 


(8) 


and this may be rewritten as 


£v _ 1 2 f* ^^t) ^0 x , x 

V Q n c Jo d0 cos 0 - cos 0 Q 

Equation ( 9 ) is strikingly similar to an integral ocurring in the 
theory of the lift distribution of a finite wing in incompressible flow. 
There , the induced angle ai is given by 




d0 * 


cos 0* - cos 0 


( 10 ) 


where 7 is the local dimensionless circulation. 


Multhopp (reference 4) has given a solution for equation (10). He 
divides the range of integration into (m^ + 1) intervals (see fig. 

below) 



with 



y 



(11) 
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and computes at the points 9 n . He assumes that 7 may he 

expanded in the form 

7 =X! c p si11 ^ 

or 


m 


7 = 


m + 


^,2 sin n 0 n sin p .0 
x n=l 


ti=l 


He then obtains the expression 

m l , 

a iv = b W 7 v " ^_i b vn 7 n 


( 12 ) 


(13) 


The prime on the summation symbol indicates that ' n = v is to he 
omitted from the summation because that special term has already been 
considered in the first term of the right-hand side (i.e., by V 7 v ). 
Reference 4 presents tables for the coefficients by V and b vn for 
mj = 7, 15, and 31- Applied to the problem at hand, m^ = 31 would 
appear to be rather small; therefore a table for m^ = 63 has been 

computed and is included in the present report (appendix A). As a 
comparison: For = 63 , A 9 = 2.8125°; for Naiman’s method with 

160 points, A 9 = 2 . 25 °. 

Utilizing this method of integration which was developed by means 
of Fourier series, an expression may be obtained for the velocity 
distribution as follows: 



(14) 


The great advantage of this method is that of simplicity: (l) The 

act ual computational procedure is very simple and ( 2 ) the derivative 


dx 


is avoided. 


The simplicity of computation is reflected in the 


fact that the time required for computing —■ at one value of 9 Q 

*0 

is approximately half that required by the method of Naiman when the 
intervals have approximately the same size. It should be noted. 
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however, that the accuracy of the method of Naiman will he greater than 
that of Multhopp in those cases where the differentiation of Ay^ hy 
Fourier expansion (equation (12)) does not give good results. 

A third method of evaluating the Poisson integral became known 
during the course of the present investigation. In a paper hy T imman 
(reference 5 )> the integral is studied in the form 

T ( 0 ) = - ~ / h(t) cot 0 ■ " — di{f ( 15 ) 

2jt Jo 2 

Timman assumes that ^(ijr) is not given analytically, hut only at 
equidistant points. An interpolation polynomial (reference 6 ) for ^(ijf) 
is employed, and these polynomials replace the function F(ijr) in a 
single interval hy a function of third order. The polynomial function 
thus introduced has a continuous first derivative,^ and it is evident 
that this continuity is essential for the attainment of good results. 

Timman has divided the period 2it into 36 intervals of equal 
length and established a computing scheme. The function crOlO is 
separated into its symmetrical and unsymmetrical parts so that 

dOlO = s + d ( 16 ) 


Then 


T (^) 


18 

= S “W 8 ! 

k=0 


18 

-2Z 

k=0 




(17) 


where the factors and 0^ 8X6 given in tabular form. In the 

d (Ayt) 

dx 


present particular case 


is antisymmetrical (equation ( 3 )) and 
. Thus the separa 

equation (l6) does not require any additional work. 


~ is symmetrical (equation (4)). Thus the separation indicated hy 
V o 


Timman *s method should give good results provided that a sufficient 
number of intervals are taken - the division of 36 intervals over a 
period of 2 it (i.e., 18 intervals over the chord of the profile) appears 
to he insufficient for an accurate representation of the function which 

dfAyt) Av 

occurs, or =— . 

7 a* 

%he polynomials used in the classical interpolation formulas 
less smooth (see reference pp. 7 and 10, figs. 1 and 2). 


are 
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The time required for computing one point hy the method of Timman 
is approximately the same as' for Naim an *s method with the same interval 
size. 


Other methods of evaluating the Poisson integral have been 
suggested. They will not "be discussed here as it is the intention of 
this section to consider only the most practical of the known methods. 
The three methods already discussed have their own particular 
advantages and have been especially developed for rapid and simple 

computation: however, all three of these methods, when — » *- or 

dx V D 

change rapidly in magn itude, become cumbersome, and require that very 
small intervals be taken over the entire range of integration because 
the scheme of equal interval size is utilized. 


EVALUATION OF POISSON INTEGRAL BY A METHOD 
EMPLOYING UNEQUAL INTERVALS 
Development of Method 


As the change in airfoil shape, or the change in velocity distri- 
bution, is given originally as a function of x it appears logical to 
retain the coordinate x in selecting the size of the different 
intervals. Hence, the Poisson integral may be Btudied in the form 


T 




( 18 ) 


which corresponds to equation (l). Conforming with its physical 
d(£y, ) 

meaning cx(x) = — — — is assumed to be a function which is finite in 

dx 

every point of its range of definition.-* 

Define 


= x n+l - *n 


(19a) 


with 


n — 0, 1, 2, 3 , .. . 


^This restriction will be dropped later; see discussion beginning 
with the first paragraph after equation ( 32 ) . 
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and 


*m < x o < 


x 

m+1 


For convenience , there is chosen (see following fig. ) 


jjm + Xm+l 
2 


(19b) 


(19c) 


0 


+ 







C 


The function cx(x) is approximated hy straight- line segments (see third 
sketch in preceding section) . Then, for x Q < x < x ja+1 , 

cr(x) = p^XqJ + — M (x _ xj 


°n + 


J n+1 


- a 


AXq 


”( x - *n) 


( 20 ) 


from which there is obtained 
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by definition. Tbe function j no ^ in the different regions of x, is 
given by different expressions as follows: 


4no 


X n+1 X ° ^ 

log e — for x n + i> x n> x c 

•*n A o 


^ x n+l ‘ *o _ 

ioge _ — for xn +1 > Xq > ^ ^ 


o - *n 


Xq 3C_ . 1 

log e — — for x o > Xn+l > Xq 


(23) 
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Introducing j no into equation (21), there results 


T (^o) = - 1 


= _ 1 
Jt 

Or, defining 


.a) + V 

ff + Vl “ a nr. 

n l + 2_, 

n Ax n (' 




J no 


Jno , 


(24) 


l + . 5>JL*a 3 

Ax„ J no 


= J 


* 


no 


(25) 


there results, finally. 


T <*°> - - i 


°n^no + | ( CT n+l " °n) 


J no 


(26) 


Since x n+1 = x^ + Ax^, the functions j no and j nQ * may he 
written as 


j * „ 1 + *o - j 
J n° Ax n 


Jno = log e ^- + ~ n ^ n Xo ) for *n > -o 


= l0g e^- + for x n + Axn > x Q > x Q 

Ax n 


= l°ge 1 + 


(27a) 


Y *u - V f ° r ^ > "n + ^ 

and this form shows that j Q0 and j no * are functions of -S-I *0 


only. 


Ax n 
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For 


*n 


- x. 


AX r 


For x Q - x n 


= 


J no 


J no 


0 


= 0 


1 0 
d no 7 


1 * = 1 
d no 


For very large 




Ax r 


= 6, 


J no 


Jno 


* 2& 


+ • • . 


2i 


3 5 


+ » 


For very large negative 




- x„ 


**n 


with 


Xn " x Q 




= i*. 


“no 


Jno 


e * -*2 

5 26 


* „p 

26 36 *^ 


(2Tb) 


(27c) 


( 27 d) 


These functions are given in figure 1 and in table I. 


It is seen that -high absolute values, of j no and j nQ occur 

X n ~ x o 4 

near those values of — — which characterize the critical interval. 


Ax n 


Figure 1 gives an idea of the characteristic qualities of j no 

y . ^ 

and j no as functions of — ; however, the representation is not 

Ax n - - 

sufficient for picking out values for a computation. Table I gives the 

y _ 

values of j nQ and j no * for -49.5 < — — < 49.5. This table might 


\f x Q = — * XnH ~l (equation ( 19 c)) the critical interval is 


given by -0. 5 < Xn ~ — < 0. 5- 

Axn 
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be used for rough computation and for getting acquainted with the 
method. In general, it is advisable to use those tables which are 
given in appendix B. 


It will prove of benefit to investigate the exactness of that 
portion of the integral which contains the singularity. Recalling 
that the function cr(x) was replaced by a straight line in every 
interval (equation (20)), there is obtained: 


f x o + lT 

* / = - 7j(<w - %) ( 28 ) 

J X °~ 2 


if x^ < x Q < x _ . Now, let an expans ion of the function cr(x) in 
the critical interval around Xq be assumed as follows: 



Then, 



Comparison of formulas (30) and (28) shows that the error in the 
critical interval is approximately given by 
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The error of evaluating the whole integral by finite differences may 
be est ima ted by using two different interval distributions and comparing 
the results for a given x Q . 


However, in addition to that error of the result produced by 
replacing the Poisson integral by a sum there exists another error. 
This sum cannot be computed exactly, but has a certain error depending 
on the accuracy of the given data for cr n (x) and the tabulated values 

„ <a/Ay t ) 

of j no and j . As the function cr(x) = — - usually has an 

nu no dx 

error of ^ = I X 10“ ^ it has proved amply satisfactory to give J nQ 

and j no * to four decimal places, the error being less than 

e 2 = 5 x lO - ^. The error of 



is smaller than its upper limit given by 

Ur 


m + m * 2 - 4 


(32) 


This formula shows that the influence of ep is stronger than the 
influence of as long as | [0 n | + X, | |cr n+ p - cr n | is smaller 

than 1 - as it is in our later examples - and the sums Xi no I 
and XVl are always larger than 1. An increase of subdivisions 

make s the sums in the upper limit of the error (32) grow, thus requiring 
a higher accuracy, especially in cr n and perhaps also in the values 

of j no and j nQ *. 


In establishing the solution of equation (18) it was assumed 
that cr(x) is finite throughout its range of definition. If it is 
desired to compute the change in shape due to a proposed change of 
velocity distribution, this restriction must be eliminated, as will 
be recognized Immediately. 
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Equation ( 5 ) may be written in the form 


d ( A yt) 

dx 


1 

jr 




Av/V 0 dx 
Vx(c - x) x - x Q 


(5*) 


Omitting the factor \J x Q \c - x^j , which does not affect the integration 

process, the integral may be reduced to the form of equation (18) by 
defining 


== = ^(x) (33) 

x) 

However, ^(x) will be infinite at x = 0 and x = c if ^ 0 

/Av\ yoJo 

and (=— ] -f 0; therefore, a special consideration of the neighborhood 
VO/c 

of x = 0 and x = c is required. • This is done by splitting the 
integral into the. following three parts: 




with and being small compared with c. The integral 

1c_€ 2 , , dx 


may be treated as was explained formerly for 
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(see equation (l8)) "because . cr;j_(x) is finite for < x < c - gg* 

For the first and third integrals, however, a new integration formula 
must he developed.. By introducing 


(J. = c - x 


and cr^x) = a- L |i(x) = a 1 (u) 


there is obtained 


ct-, (x)- 


dx 


*C -6, 




2 *f\ < 3-4 

<*1 ( m -) 

10 - ^o 


( 35 ) 


Hence, the method used for the first integral will also apply to the 


third. In most cases 


will he zero and there will he no need 


for a special evaluation in the neighborhood of x = c. 


The- integral 




Av/Vp dx 

\jx{c - x) x “ x o 


( 36 ) 


will have an important influence on the result of equation (3*0 only 
if x Q is near to eq. First the general formula will he given and 

then a simplification will he discussed for x Q » e-^. 

The integral ( 36 ) will he solved assuming that 

f; = a o + & i(f) + s 2 (f) for 0 <• x < € i 

Only the final formula of this procedure is given here; the 
details of the solution will he found in appendix C. 


( 37 ) 
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with M 0 given In figure 2, and 


with 


with 



( 39 ) 


The coefficients a Q , a^, and a^ may he determined first, as 
they do not depend upon the particular value of x Q , and then Fj^x 0 ^ 
may he computed. The term depending on a^ and a 2 will exert an 
influence only for small values of x 0 /c. After a brief training the 
computor should he able to decide rather accurately when the formula 



is sufficient and when the more exact expression (equation (38)) is 
required (also see fig. 2). 
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Organization of Computational Procedure for Unequal Intervals ; 

Transition from One Size of Interval to Another 

A thorough, underst andin g of the method is best achieved by 
following through a rather simple example; in addition various short 
cuts to the method will be demonstrated. 

Assume a function cr(x) of the type shown in the following 
figure. 



It appears reasonable to take rauher small intervals for small values 
of x because of the form of the curve cr(x); therefore, the following 
arrangement of interval sizes is arbitrarily selected: 


Ax = 0.002 for 0 < x < 0.030 

Ax = 0.006 for 0.030 < x < O.O 96 

Compute t ( 0.009) with the help of equation ( 26 ). Note that the 
critical interval, extends from 0.008 to 0.010. Table II (a) gives 
the values of x/c, °n+l " V and J no* for ^ range 

with Ax = 0.002. At x = 0.030 the interval changes to Ax = 0.006 
and the same functions are given for the range with this size of 
interval in table II (b). Naturally the range above the broken line 
in table II (b) is not utilized in the computation since this portion 
has been considered in table II (a) « 
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x - x p 

Note that progresses in table II (b) in the same manner 

as in table II (a); this is due to the special choice of Ax. If 

Ax = 0.006 were used starting with x = 0 the critical interval for 

x Q = 0.009 would extend from 0.00 6 to 0.012. Hence, for Ax = 0.006, 

* 

jno = 0 and j no = 1 is to be found at x/c = 0.006. 

For rapid computation it is best to have j and j n * as 
functions of ^ ^ 0 on a paper strip and to place this strip adjacent 

X ""X 

to the columns headed by cr n and cr n+ j_ - c n . If n ■ — 2. progresses 

as indicated in table I, the correct location of j no = 0 and j no * = 1 

at the beginning of the critical interval fixes the placement of the 
strip. 

In the example just treated, the transition from one size of 
interval to another is very easy because x Q lies at the midpoint of 
an interval of the size 0.006 as well as of the size 0.002, if starting 
with x = 0. 

If Ax had been chosen 0.004, such a desirable arrangement would 
not have resulted because x Q = 0.009 would not be located at the mid- 
point of an interval of this size (starting with such intervals 
at x = 0). 


As a second example compute the value of t at x Q = 0.015. 

Xn - x Q x n - x Q 

Again, — — and — — will progress as in table I. The values 

Ax Ax 

j no = 0 and j QO * = 1 will be placed opposite x/c = 0.014 for the 

region with Ax - 0.002 and opposite x/c_= 0.012 for the region 

with Ax = 0.006. As long as Ax = 3Ax, Ax = 3Ax, and so forth and 
if x Q is chosen so as to be at the midpoint of the largest size of 
interval, the computation may be accomplished by shifting the strip 
with j no and j no * corresponding to table I. 


But suppose that the interval sizes are so arranged and it is 
desired to compute a point where x Q does not lie at the midpoint of 

the largest size of interval; for example, x Q = 0.013. The value 

x Q = 0.013 lies at the midpoint of an interval with Ax = 0.002; hence. 
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for the range 0 < x < 0 . 030 , j no and j nQ * may be taken directly 

from table I. However, at x/c = 0.030, intervals of the size 
Ax = 0.006 commence and there is obtained 

*n ~ X Q = 0.030 - 0.013 = 2< 833 
^ 0.006 


The value of ^2 ^2. progresses by 1, that is, 2.833; 3.833, 

Ax 

4.833, .... Thus the functions j no and j n0 * are needed for 

values of — — which are not given in table I. One might think of 

Ax 

taking them out of an enlarged diagram (see fig. l); however, it iB 
much more convenient to take them out of an extended table, which is 
conveniently arranged for ‘'advancing by 1." Such tables are given in 
appendix B. 

The example presented by the figure at the beginning of this 
section suggested starting at x = 0 with the smallest intervals. 
However, other examples may suggest another distribution of intervals. 
The smallest size of intervals may lie at any part of 0 < x < c. 

There are no restrictions in the arrangement of intervals. (See, e.g. , 
discussion following equation (43).) 


Accuracy of Method, Examined by Means of 
an Analytical Example 

The accuracy of the result depends directly upon the size of the 
interval taken and the reliability of the data comprising the func- 
tion a(x). Because the function c(x) will be replaced by a broken line, 
a glance at the curve will quickly Suggest an arrangement of intervals. 

In addition, the error in the critical interval may be used as a first 
test of the choice of intervals . . 

As a test of the quality of this new method, involving unequal 

a(Ay t ) 

intervals, a function a(x) = — has been treated which allows 

dx 

the analytical computation of t (x) = . 

* o 
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The function cr(x) is given analytically as 


0 < 2A 


d ( A yt) 


= Bx(2A - x) 


2A = x c^ • — ■ = -D^ - x)(x - 2A) 


C]_ < x "h c 


d fot) _ 


The following arbitrary values have been selected: 


c-l = 0.35 


c = 1.0 


D = Some multiple of B so that 


'! d(Ay t ) 


dx = 0 


The functions Ay + and — — 42. are given in figures 3(a) and 3(b), 

dx 

respectively. 
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At 


The analytical computation of =— for figure 3(h) is given in 

’ o 

figures 4(a) and 4(h). The arrangement of the unequal division for the 
numerical computation of is indicated in figure 4(a). ^ 

* n 


^ It vas desirable to obtain the value 
the first interval has heen placed so that 

d fot) „ 


function 


dx 


for -0.001 < x < 0. 


of QL at x n = 0; hence, 
v o 

-0.001 < Xq < 0.001 and the 

d 

Since the function — 

dx 



is replaced in every interval hy a straight line, the error might he 
expected to he large. However, g = 0 at x Q = 0 will aid in 
preventing the error from "being too large. 



A more exact solution would he obtained hy putting 


g = 0 - ^ <x <0 

g = ~x 0 < x < — 
Ax £ 

2 


For the interval 



Ax - 
— < x 


<Ax 


equation ( 26 ) would yield 


. 1 r to/2 g _£l_ 

*j- Ax/2 x - x ' 


1 


(si - 0 ) *)no + 0 x >3no 



and no error is introduced. For x Q ^ 0 there is a very small error 
which may he avoided hy respecting the change of size of the interval 
near x = 0. 
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Also given in figure b(a) are points of the ^ curve determined by 

V 0 

the method of unequal intervals. Figure b(b) presents the same 
lnfonnation plotted to a larger scale. 

For comparative purposes the same problem has been treated by the 
three methods of computation discussed earlier, namely, those of Naiman, 
Multhopp, and T human . Figures 5(&) and 5(k) show the results obtained 
by the method of Naiman; obviously, the bO-point solution does not use 

dJAy+) 

a sufficiently accurate representation of the — ^ curve, while the 

dx 

80- and 160-point solutions are quite good, with the exception of the 

maximum and minimum points of the curve. In order to obtain a 

* o 

value at approximately ^ = O.O36 a solution involving 320 points would 

be required. In this respect the method of unequal intervals is more 
adaptable to special conditions without involving much new work than is 
the method of Naiman. 


The results obtained by Multhopp’s method are given in figures 6(a) 
and 6(b). The 31-point solution (in Multhopp’s somewhat odd manner of 
designation) corresponds to Ad = 5*625°; the 63-point solution, to 
A0 = 2.8125°. The computation is very simple and the results of the 
method with 63 points are comparable with that of Naiman with 80 points, 
with the exception of those near the region 0 < x < 0.01 (this iB 

shown most clearly in fig. 6(b)).' The very Bteep peak of — 

dx 

at x/c =0.02 requires rather high harmonics for the representation 
of Ay±i consequently, good accuracy in the region near the origin may 
not be expected. This is substantiated by the fact that for the 63-point 
method the highest effective harmonic would have three waves in the 
region 0 < x < 0.0b; obviously a sufficient degree of accuracy in the 
differentiation process cannot be obtained. 


As mentioned earlier, Timman’s method might be expected to give 
good results if the size of interval is properly chosen. Inasmuch as 

= 3§0 = 1q0 


only a table for A 0 = 


computation for — 
• o 


= xu- was available, the result of the 
cannot be expected to be good, as is evidenced by 


observing figure 7. Ike result obtained is comparable with that of 
Multhopp’s 15-point and Naiman *s bO-point solutions. 


An excellent method of examining the accuracy of these methods 
still further is simply that of solving the inverse problem. From the 
av „ _ d (Ay t ) 


curves of just discussed, values for 

Vo ' dx 


have been computed 
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and are presented in figure 8. The method of unequal intervals gives 
good results, indicating that the arrangement of Intervals chosen was 
as good for the inverse problem as for the direct problem. It is 
apparent that Naiman’s method requires even smaller divisions than 
l60 points in order to avoid inaccuracies near the point x/c = 0 . 04 . 

The reader may wonder that the inverse problem is not given by 
Multhopp 1 s method. It must be recalled that Multhopp*s method of 
solving the direct problem does not involve the differentiation of Ay^; 

that is, it is particularly fit for this problem and presents, on the 
other hand, no analogy for the inverse problem: 


= A l' 2 * cot 0.9 . - i r* mis ds 

dx 2« J 0 V Q 2 it Jo V Q cos 6 - cos 9 0 

Because more extended tables for Timman’s method are not available, 
and the results obtained from the 36-poinh method for which tables exist 
are very poor, no further examples of the application of this method 
will be given. 

COMPARISONS OF METHODS OF NAIMAN AND MULTHOPP WITH METHOD EMPLOYING 
UNEQUAL INTERVALS BASED ON ACTUAL EXAMPLES OF CHANGES IN 

AIRFOIL SHAPE 


The method of unequal intervals has shown good qualities when 
applied to a problem where the function cr(x) is known analytically. 
However, as mentioned earlier, this function is not usually known in 
analytic form. This section, therefore, will compare the three 
principle methods, those of Naiman, Multhopp, and unequal intervals, 
on the basis of actual design problems, solving the direct problem 

for ~ and using these results to solve the inverse problem (excluding 
*0 

Multhopp for the inverse problem). 


Figure 9 (a) shows the Ay^ relations for examples I and II and 
d(Ay.) ri/A " ' 

figure 9 (b), the ■ — relations. Note that the slope of 

for example II is more than twice that of example I near x/c = 0 . 


11 ana 

PKl 
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The direct problem for example I by Naiman’s method is given in 
figure 10. The l60-point solution does not show any appreciable 


deviation from the 80 -point solutions at the region of 


however, near the origin, at 


/Av;\ 

Vol 


'max 

the influence of the smaller-sized 


'min 


intervals (8O:A0 = 4.5°; l6O:A0 =2.25 ) is quite pronounced. 

The solution by Multhopp ’s method is given in figure 11; ^ 31 points 
around the half circle are not sufficient for a solution comparable with 
Naiman’s 80-point solution, and even a solution based on 63 points does 
not offer much improvement. The results are poor, as might be expected, 
in the region very near the origin (see preceding section). 

Figure 12 presents the results obtained by the method of uneq ual 
intervals, compared with results obtained by Naiman’s 80- and 160-point 
solutions. The method of unequal intervals gives results corresponding 
to those established by Naiman’s l60-point solution. The subdivision 
used is shown in the figure. 

As before, the inverse problem was solved, and is given in fig- 
ure 13. In each case the computed curve of was the one used in 

d (Art) 

dx 

results, thus proving that the chosen number of divisions was sufficient 
in Naiman’s method and in the method employing unequal intervals. 


obtaining the values for the 


curve. Both methods give good 


The value of 


computed at x/c = O.I 71 is of some interest. 


d fot) 

dx 

This point was computed by the method of unequal intervals in two 
different ways: First, the arrangement of intervals shown in figure 13 

was utilized to compute the lower point. Then a new arrangement of 
intervals (Ax = 0.018 for 0 < x < O. 36 ) was set up and the same point 
computed. The idea was to determine the inaccuracies that would result. 
One might predict that, since the point x/c = 0.171 lies at a con-' 
siderable distance from the region of rapid c han ges in errors of 

Vo 

only small magnitude would be introduced; this is fairly well sub- 
stantiated by the results shown in the figure because the error thus 
introduced is approximately that of the deviation of Naiman’s l60-point 
solution. 


^Recall 
of Ay t . 


that this method does not involve the differentiation 
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Nov, turning our attention to example II, which, it will "be 


recalled, has a slope of 


L (Art) 

dx 


of approximately twice that of 


example I, the results given in figures 14 to 17 are obtained. 


For the direct problem Naiman’s method of 160 points and the 
method employing unequal intervals give results which are in good 
agreement. For the inverse problem (fig. 17) it is apparent that the 


method using unequal intervals is superior see the deviation at 


L (Ay t ) 


dx 


Multhopp *s method gives 


given by the method of Naiman 

— 1 max 

a rather good result (fig. 15 ), which may be attained when the Fourier 
representation of the Ay curve is adequate. 


The two examples thus far presented are favorable for Naiman’s 
method because the steep slopes of o(x) occur near x = 0 whifre the 
points Naiman uses are close together. However, going to still steeper 
slopes near x = 0 would require a rapidly increasing number of points. 
The new method offers another possibility here. Assume that in that 
critical region x fc < x < x k+1 (x k may be 0) cr(x) may be represented 

by ff(x) = ^ a n x n . Then the integral 



o(x) 
x - x Q 


dx 


may be split into three 


integrals 





(41) 


The first and third of these integrals may be solved in the usual 
manner using the functions j nQ and J no ** The second integral will 


be solved analytically. 
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This simple form, due to the use of the coordinate x in the 
Poisson integral, allows a rapid integration, because the integral 



can he solved by recurrence as follows: 

x n _ x n 

Vo = + *okn-l,o for 11 ^ 1 (^2) 

} n 7 


with 


k 


o, o 


= i 


no 



(* 3 ) 


Thus even very steep slopes cause no difficulties. 


As already mentioned, examples I and II correspond well to the 
qualities demanded by Kaiman’s method insofar as the rather steep 
slopes occur in those portions where the points 0 n are close together. 
If those steep slopes should occur in other portions of the chord, how- 
ever, a very great number of points in the Kalman method would be needed 
in order to represent d(x) adequately, and to get reliable results. 

In such a case the method using unequal intervals shows its advantage 
by allowing a free subdivision of the chord. 


A third example will serve to illustrate this. Figure 18 shows a 


function cr(x) 


-(Ayt) 

dx 


The essential values of the function lie in 


a part of the chord where even Kaiman’s method with 160 points is not 
sufficient to represent the function accurately. This is forcibly 

shown by the two curves of If the function cr(x) is modified 

V ° A 

(dotted line) so as to eliminate the high peak, then the ^ curve by 

v o 

unequal intervals can be made to agree with the original ~ iby Kalman 's 

*0 

160 -point solution, thus definitely proving that, in this example, 
Kalman’s method with 160 points is insufficient. 


Table III indicates the computation for the point x Q = 0.0 65 
by unequal intervals. 
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CONCLUDING DISCUSSION 


The new method of evaluating the Poisson integral developed herein 
is to be reco mme nded for all those functions cr(x), where steep slopes 
in small portions of the region to be integrated exist. In these 
portions a very small size of interval may be chosen without requiring 
that this same size of interval be used throughout the region of 
integration. In this manner, the work required for computation may be 
maintained at a reasonable level even for the most complicated problems. 

The anal ytical treatment of special parts of the integral is 
possible (evaluating the remainder by the new method; see preceding 
section). In those problems where a transition to very small intervals 
in part of the integration range would require the determination of a 
great many values of cr n , this idea might be used to advantage. 

It should be noted that the smoothness of the function a(x) and 
its accurate representation by single points is essential for good 
results. If, for example, single points cr n are simply taken from a 
curve for x n very close to one another it may be compulsory to check 
these values by a table of differences. 
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7.3835 

.7760 

♦ B677 

.1317 

.0774 

.0907 

43 

25 

.0458 

,u48l 

, OdO-a 

,06/3 

,0584 

.0855 

.1180 

. 180 7 

.3555 

.8613 

y.oooo 

7 * v £74 

.7509 

• CO.L-Q 

• ii*W 

•0773 

41 

as 

.0340 

.0366 

.0303 

,0427 

,0490 

.0601 

.0774 

,1074 

.1668 

.3078 

.8146 

7.0388 

3.7703 

.7318 

.2373 

.1895 

39 

87 

.0858 

.0270 

.0397 

.0337 

.0371 

.0438 

,0639 

.0704 

,0969 

.1661 

.8906 

.7781 

0,7364 

6,0108 

.7907 

.2966 

37 

£9 

.0315 

,0222 

.0233 

.0265 

.0285 

,0330 

.0394 

.0492 

•0660 

.0985 

.1467 

.2777 

.7496 

5.5335 

6,4941 

.7173 

53 

31 

.0174 

.0180 

.0109 

.0204 

,0225 

.0256 

,0298 

.0360 

.0455 

.0608 

.0873 

.1402 

.2679 

.7313 

0,3297 

6. 30 £3 

33 

33 

.0143 

.0147 

.0164 

.0156 

,0180 

.0201 

,0831 

.0873 

.0334 

.0487 

,0076 

,0037 

.1363 

.*ai« 

.7828 

6,5019 

51 

53 

.0110 

.0122 

.0137 

.0166 

.0146 

.0158 

.0103 

.oaiB 

,0853 

.0313 

.0404 

.0551 

.0009 

.IMS 

.2501 

.6149 

89 

37 

.0099 

.0101 

.0106 

.0111 

.OlfiO 

.0131. 

.0147 

.0168 

.0197 

.0238 

,0297 

.0587 

.0333 

.0790 

.1304 

,7360 

27 

39 

.0083 

.0066 

.0068 

.0093 

.0099 

.0106 

.0U9 

.0130 

.0166 

.0180 

.0223 

.0284 

,0374 

.0620 

.0778 

►1296 

23 

A 1 

.rmffj 

.nn*l 

.(VITA 

JY1T7 



.nrxiA 

.min 

-M an 

*nii.A 

-M 

• run n 

.njwa 


• fWlS 


aa 

43 

.0069 

.0060 

.0052 

.0056 

.0059 

.007* 

.0081 

.0090 

.0102 

.0117 

.0138 

.0167 

.0207 

• 0855 

.0368 

,0507 

81 

43 

,0050 

.0050 

.0068 

.0054 

.0087 

.0062 

.0067 

.0074 

.0003 

,0095 

,0110 

,0131 

,0100 

. 0*01 

.0861 

,0354 

19 

47 

•004B 

.0048 

,0044 

,0045 

.0046 

.0031 

.0056 

.0051 

.0056 

.0077 

.0089 

.0104 

.0126 

,0154 

.0196 

.0287 

IT 

49 

.ora 

.0035 

.0030 

.0038 

.0040 

.0042 

.0046 

.0060 

,oow 

.0068 

.0071 

.0083 

.0009 

.0130 

.0149 

.0191 

13 

81 

.0089 

.0089 

.0030 

.0031 

,0033 

.0036 

.0037 

.0041 

.0043 

.0050 

,0067 

.0055 

.0078 

.0093 

.0115 

.0144 

IB 

BS 

.0013 

.0084 

.0024 

,00 EC 

.0085 

,0088 

.0030 

.0053 

.0035 

.0040 

,0046 

.0052 

.0051 

* .0072 

.0060 

,0100 

11 

BB 

.0018 

.0019 

.0019 

.0080 

.0081 

.0082 

.0024 

.0035 

.0030 

.0031 

.0035 

♦0040 

.0047 

.0066 

,0055 

.0088 

9 

57 

.0014 

.0014 

.0015 

.0015 

.0018 

.0017 

.0018 

.0019 

.0021 

.0023 

.0025 

.0030 

.0030 

.0041 

.0049 

.0069 

7 

59 

.0010 

.0010 

.0010 

.0010 

.0011 

.0012 

.0012 

.0013 

.0013 

.0015 

,0016 

,0081 

.0024 

.0020 

.0033 

.0040 

5 

51 

,0005 

,0006 

.0005 

.0005 

.0005 

,0007 

.0007 

.0008 

.0009 

.0010 

.0011 

.0012 

.0014 

.0015 

.0010 

.0083 

3 

63 

.0008 

.0008 

.0008 

.0008 

.0002 

.0008 

.0008 

.0005 

.0003 

.0003 

.0004 

.0004 

.0000 

.0006 

.0005 

.0008 

1 

■ 

58 

60 

38 

35 

64 

9£ 

50 

48 

.46 

H 

42 

40 

38 

35 

34 

3ft 

V 


TS+iS m vdvh 
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APPENDIX B 


VALUES OF AND j * AS FUNCTIONS OF - 

D.U uu /\y 

The values of the functions j no and j no * are presented as indi- 
cated in the following table. The values are tabulated in a form selected 
to minimize the necessity for interpolation except for the region contain- 
ing the singularities of the functions j no and j^*. For ease in 
computation, tables B— I to B— VIII , inclusive, are arranged so that the 

vertical increment of X - ^° is unity. Table B— IX gives additional values 
for the region containing the singularities of the functions j no and 
Jno** 


TABLE NO. 

BANCE OF z ^ r “ s ° 

TTSjnpwmnqi QP x V~ x o 

Ax: 

Ax: 

B— I 

-189 to -90 

1.0 

B— II 

-89.5 to -40.0 

• 5 

B— III 

-39*9 to -20.0 

.1 

B-rv 

-19.99 to 0 

.01 

B-V 

0 to 19.99 

.01 

B— VI 

20.0 to 39.9 

.1 

B— VII 

40.0 to 89.5 

•5 

B-VIII 

90 to 189 

1.0 

B— IX 

-1.000 to 0.000 

0.001 




UJ 

Os 


TABLE B— I.— VALUES OF ^ AND j no * USED IE EVALUATING EQUATION (26) 

- 189 o - 90 

Ac 



-i&C.O 

-1 


-JiSt.O 

-131.0 

-1*1.0 

-13X.0 

-121.0 

-ixt.o 

— 1 QX .0 

-SS.O 

2 

Jno 

W 

Jno 

duo* 

dno 


Am 

laa* 

Jno 

■W 

mm 

Jno* 

dno 

Jno* 

duo 

3 aD* 


Jno* 

Jno 

Jno* 

9 

-0.0053 

- 0.0027 

-0.0056 

-0.0028 

-0.0059 


-0.0063 

-0.0032 





IPffiBEl 





-0,00*6 

- 0.0102 

- 0.0051 

8 

-.0053 

-.0094 

—.005*1 

-.0027 

-.0056 

-.0028 

-.0060 

3 2 S 

R§ 

Rw 

-.0032 

-.C 068 

-.003* 

-.0073 

-.0036 

-.0073 

b ffiy 

-.0065 

-.00*3 

fm 

-.00*7 

-.0103 

-. 010 * 

-.0105 

-,0051 

l 


-.0057 

-.0057 

B ES 

-.0060 

-.0060 

■ SJ 
5 22 

-.0032 

-.0032 

-.0C68 

-.0069 

-.003* 

-.003* 

-.0073 

-.0074 

-.C 037 

-.0037 

-.0079 

-.oofio 

Eft* 

— ,0066 
-.0067 

-. 00*3 

—.0043 

-.00** 

-.00*7 

-.00*7 

-.0052 

-.0053 

? 

-.00511 


-.0057 

-*0056 


-.0061 

! 2 E 

-.CO63 

-.0032 

-.0069 

-.0039 

-.007* 

-.0037 

— .cofio 

P “ 

-.0007 

-.0060 


::S! 

-.0106 

-.0053 

4 

-.005* 

S&2I 

P*jf*Vy 

-.0061 

1 vES 

-.0065 

-.0033 

-.0070 

-.0033 

-.0075 

-.0037 

-.003fl 

-.0081 

—.0041 

-.00** 

Hwf 

-.0107 

-.0106 

-.0053 

-.0034 

3 

-.0055 

n 

-,0058 




-.0066 

-.0033 

-.0070 

-.0033 

-.0075 

-.co£e 

-.00*1 

-.0069 

-.00*5 


-.00*9 

2 

-.0095 

-.0058 

-.0029 


s ss 

-.0066 

—0033 

-.0071 

-.0039 

-.0076 

-.0038 

-.cc 8 a 

-.00*1 

-.0090 

-.00*9 

-.0099 

-.00*9 

D( 4 rd 

-.0055 

1 

-.0059 

— *OCE 0 

-.0059 

—.0029 


• «^ig 

-.0066 

-.0033 

-.0071 

-.0036 

-.0077 

-.0038 

-.0083 

-.ooeS 

-.ode 

-.0090 

-.00*9 

-, 0 K >0 

-,oo?o 

-.0110 

-.0055 

0 

-.0056 

-.0028 




■ uffl 

-.0067 

-0033 

-.0072 

-.0036 

-.0077 

-.0039 

-.oo*a 

-.0091 

-.00*6 

'-,0101 

-.0050 

-.0112 

-.0056 


table; B-n,- VALUES of AND j no * USED in evaluating equation ( 26) 

- 89.5 = - to.o 

Ac 



-$L 5 

- 66 C. 0 

- 71.5 

-71.0 

-6CC.3 

-61.0 

- 51.5 

- 51,0 

. -ta.s 

-*i.O 



Jao* 

Jno 

duo* 

dnn 

d»0* 

Jno 

Jno* 

dnO 

Jno* 

drn 

dl>0* 

dao 

Jno* 

dno 

Jno* 

dnn 

Jno* 

duo 

Jno* 

9 ' 

-0.0112 

-0,0056 

-0.0113 

-0.0057 

- 0.0127 

t-0.006* 

-0.0127 

- 0 .C 06 * 

-0.01*3 

-0.007J 

-0.01*6 

— 0.0073 

-0.0169 

-0.0065 

- 0.0171 

- 0.0065 

-o.oao* 

-0.0102 

U).0206 

- 0.0103 

8 

-,ou* 

-.0057 

-.0114 

-.0096 

-.0120 

-.006* 

-.0129 

— 006* 

-.01*7 

-. 007 * 

-.01*8 

-.0074 

-.0172 

-.0066 

-.017* 

-.0087 

-.ceo6 

-.010* 

-.0211 

-.0105 

7 

-.0115 

-.OO5O 

-.0116 

-.0050 

-.0130 

—0063 

-.0131 

-0069 

-.01*9 

-.0075 

-.0150 

-.0075 

-.0175 

-.0068 

-.0177 

-.008 9 

-.0213 

-.0107 

-.0215 

-.0108 

6 

-.0116 

-.0058 

-.0117 

-.0099 

-.0132 

— .0066 

-.0132 

— 0067 

-.0159 

-.0076 

-.0153 

-.0076 

—.0179 

-.0090 

-.0100 

-.0090 

-.CC 17 

-.0109 

-.0220 

-.0110 

? 

-.0118 

-.0059 

-.0110 

-.0059 

-.0133 

—OO67 

—013* 

-.0067 

-.015* 

-.0077 

-.0155 

-.0070 

-.0102 

-.0092 

-.010* 

-.0093 

-.0222 

-.0112 

-.0225 

-.0112 

* 

—.0119 

-.0059 

-.0120 

-.0060 

-0135 

-.0067 

-.0136 

-0060 

-.0156 

-.0078 

-.0157 

-.ooj9 

-.0189 

-.0093 

-.0167 

-.0094 

-. 0 C 27 

-.on* 

-.0230 

-.one 

3 

-.0120 

-.0060 

-.0121 

-.0061 

-0037 

-.0066 

-.0138 

-OO69 

-.0159 

-.0079 

-.0160 

-.ooso 

-.0189 

-.0094 

-.0190 

-.0095 

-.0233 

-.0117 

-.0235 

-.0118 

2 

-.0122 

-.0061 

-.0123 

-.0061 

-0039 

—0070 

-,oi 4 o 

—0070 

-.0161 

-.0081 

-.0163 

-.0081 

-.0192 

-.0096 

-.019* 

-.0097 

-.0238 

-.0120 

—.02*1 

-.0121 

1 

-.0123 

-.006s 

-.0124 

-,006a 

-.0141 

-.0071 

-.01*2 

— 0071 

-.016* 

-.0082 

-.0165 

-.0063 

-.0196 

-.0099 

-.0196 

-.0100 

-.0244 

— 0122 

~.CG *7 

-.012* 

0 

-.0125 

-,ooSe 

-.0126 

-.0063 

-.01*3 

-.0071 

-.01** 

— 0072 

-.0167 

-.006* 

-.0168 

-.0004 

—.0200 

-.0101 

-.OCOG 

-.oioe 

-.0290 

-.0025 

-.0253 

-.0127 
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TABLE B— HI.— VALUES OF J no AND j no * USED IN EVALUATING EQUATION (2 6 ) 

- 39.9 = 3c5a $ -£0*0 


\I^ 


9 


B 


7 


6 


5 


4 






L 



Xrr-oco 

Ax 

Jno 

Jno’ 

4do 

Jno’ 

4oo 

to’ 

4no 

Jno* 

^DO 

^no* 

^ no 

to’ 


4oo* 

^ no 

to' 

4 bo 

tn* 

*^no 

to* 

-39-X 

33.0294 

-0,0127 


g STS 

g H5 

-0.0128 

-0.0256 


-0.0256 



-0.0129 

-0.0258 

-0.0130 

-0.0258 

^t.0130 

-0.0259 

-0.0130 

-0.0260 

-0.0131 

-38. 

-.0260 

-.0131 



||| 

-.0132 

-.0262 


-.0263 


-.0264 

-.0133 


-.0133 

-.026? 

-.0133 

-.0266 

-.0134 

-.0267 

-.0134 

-37- 

-.< 26 ^ 

-.0134 

-.0268 

||S • 

E ® 

-.0135 

-.0270 

-.013? 

-.0270 

-.0136 

-.0271 

-.0137 


-.0137 

-•0273 

-.0280 

-.0137 

-.0273 

-.0281 

-.0137 

-.0274 

-.0138 

-36. 

-.0275 

-.0283 

-.0136 

-.0276 

E • TP. 

|pf HE 

-.0139 

IjCi&jrri 

-.0139 

—.0278 

-.0139 

-.0279 

—.0287 

-.0140 


-.0140 

-,oi4i 

-.0141 

-.0282 

-.0141 

:£ 

-.0143 

-.0283 

ke 


-.0142 

-.0143 

-.0286 

-.0143 

-.0144 

-.0145 

-.0268 

-.0145 

-.0269 

-.0145 

-.0290 

-.0146 

-.0291 

-.0146 

-.0292 

1 5 SI 

-*0292 

->0147 

-.0293 

-.0147 

-.0294 

-.0148 


-.0148 


-.0149 

-.<*97 

-.0306 

-.0149 

-.0296 

-.0150 

jEIpl 

-.0150 

-33. 

-.0299 

“,0150 

- 0300 


-.0301 

-.01?! 

-.0302 

-.0152 

-.0353 

-0153 

-.0153 

-0305 

-.OI54 

-.0154 

-.0307 

-.0154 

-.0155 

-32. 

-.0909 

-.0136 

-.0310 

-.0136 

-.0311 

-.0156 

-.0312 

-.0157 

-.0313 

-0157 

-.0314 

-.0158 

-0315 

-.01J8 

-.031? 

-.0158 

-.0316 

-.015? 

-.0317 

-.0328 

-.0160 

—31. 

-.0319 

-.0160 

-.0320 

-.0161 

-.0321 

-.0161 

-.0322 

-.0162 

-.0323 

-.0162 

-.0324 

-.0163 

-.0168 

-.0325 

-.0163 

-.0326 

-.0164 

-.0337 

-.0338 

-.01 6^ 

-.0165 

-30. 

-.0329 

-.0169 

-.0330 

-.0341 

-.0166 

-.0331 

-.0342 

-.0166 


-.0167 

-■0333 

-.0167 

-.0334 

-.0346 

-.0336 


-.0337 

—.0348 

-.0169 

-.0170 

mi'i 

-.0171 

-e?. 

-.0340 

-.0171 

-.0172 

-.0172 

-.0173 

-0345 

“-0173 

-.0180 

-.0174 

-.0180 

-.0347 

-.0360 

BTiiTlI 

■■wCi 

WTnKil 

-.0175 

-.0182 

-.0350 

-.0176 

-.0182 

E F ; * 1 

-.0176 

—.0183 

-ee. 

-.0352 

-.0177 

-.01 64 

-.0353 

-.0176 

-.0164 

-.035? 

-.03® 

-.0178 

-.0185 

-.0356 

-.0179 

-.0186 

-.0357 

-0358 

-.0361 

-.0362 

-.0364 

-C7- 

-.0365 

-.0366 

-.0369 

—.0270 

-.0186 

—.0372 

-.0187 

B&f 


-.0375 


-.0376 

-.0189 

-.0377 

-r0190 

-e6. 

-.037? 

-.0394 

-.0191 

-.0380 

-.0191 

-.038a 

-.0192 

-.0383 

-.0193 

— .o5s^ 

-.0193 

—.0386 

-.0194 

-.0195 


-.0196 

-.0391 

-.0196 

-.0392 

-.0197 


-.0196 

-.0395 

-.0199 

-.0397 

-.0413 

-.0200 

-.0398 

-.0200 

— .o4oo 

-.0201 

—.0402 

—.0202 

ERSC 

-.0203 

JkIC 


-.0407 

-.0205 

-.04o6 

-,020b 

-.0410 

-.oeo6 

-.0*12 

-.0207 

-.0208 

-.0415 


-.0417 

WKlfll 

-.04 18 

-.0211 





-.0424 

-.0213 

-.0426 

-.0214 

-23. 

-.0447 

-.0215 

-.0439 

-.oai6 

-.0431 

-.0217 

-.0433 

-.0433 

-.0219 

-.0437 

-.0219 

-.0439 


SKJjJP 


-.0443 

-.0223 

-.0444 

-.0223 

-22. 

-.0447 


—.0449 


-.0451 


-.0453 

-.0228 

-.0455 


-.0457 

-.0230 

-.0459 


SBt!y 

-.0232 

-.0463 


mm 

-.023? 

-CL. 

-.0467 

-.0236 

-.024T 

-.0470 

-.0237 

-.0248 

-.0472 

-.0238 

-.0474 

-.0239 

-.0476 

-.0240 

-.0479 

-.0241 

-.0481 

mm 


-.0243 

-.0486 

-.0246 

-eo.x 

-.0490 

-.0493 

-.0495 



-.0251 

-.0500 

-.0252 

-.0503 

-.0253 

-.0505 

B m 

-.0256 

-.0510 

— 02?7, 

-0513 

—.0259 




i 
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TABLE B— IV.— VALUES OF AMD USED IX EVALUATING EQUATION (26) 


(jj 

CD 


(a) - 19 .X? S < _o <xx vhfire 99 ^ xx > 90 

Ax ' ' ' 



(l) -19.22 < ^ y C ° $ -O.XX where 89 £ 22 £ 80 

Ax 



J j j j 

































IABIE B— IV.— COHTINUED 



NACA IN 
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TABLE B— IT.— CQNTJJNUJSD 


(e) -19.XX = < —0,22 where 59 £ 22 £ 50 

As 



TABUS B-IT.- C0HTINUEU 


(f) -19.XX = < — o.SX where .49 £ 32 > 40 



idl VOVM 
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TABLE B— IT.— CQNTJJNUJSD 


(e) -19.XX = < —0,22 where 59 £ 22 £ 50 

As 



TABUS B-IT.- C0HTINUEU 


(f) -19.XX = < — o.SX where .49 £ 32 > 40 



idl VOVM 

















































































































TABLE B-V.- VALUES OF Jn 0 AHD j no * USED IN EVALUATING EQUATION (2 6 ) 
( a ) O.ZK = = 19. XX where 00 = XX = 09 


H-— -I v r, Jr V | ^» |*»* ~*» | W [w jy I *.•' ^ 

°‘SS °-**T 3 *oMo 0.6* t« t.»UT o.€rrr 2-7169 ' 


•0*3 .0319 .06*3 

.0606 .0300 r060£ 

•vp .oaBa .orn 
.Wi .0668 .05*0 

•P3U ,«* ^513 




TABLE B-V. - CONTINUED 

(t>) O.XX ~ - n ~~ — ° a 19.ZX where 10 = 3 CC = 19 


B*l» 0 . 7 lfl 9 
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r— 

5 - 


(c) o.ZK < 19 - 2 Z where 20 < XX 5 29 



TABLE B-V. — COHTIHUED 

(d) O.XX < Xn L Z ° - = 19 .XI where 30 ^ XX ^ 39 
At 
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iff A «e5K«WJisps 


TABLES B-Y. - COHTIMJEi) 


(©) 0.52 = Zg ~ X ° - 19 .52 where 40 < 22 < k-9 

Ay 
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TABLE B-7. - CONTINUED 


0.2Z « — 5 19.21 -whore 80 ^ XX ^ 89 



^rta ^oo* 
CUilM 00515 

: Mm 

•W a44$ 
J336 aifij 
afipi .0917 
U55a .077* 
•137a .0670 

.1406 .0591 

.1076 .0583 

:S§? S 

.oeu ,o*ca 
.OWE .0371 
■owe ,o}*£ 
-.063* .pjsa 
•«u .ojo* 
.owe .aefl6 

.«*7 *0B7X 


0.7900 0.3*3® 

.*>50 ,aas 

.JMS .1*37 
JtJSO .1115 
nftftp .0911 


aace .0589 

3g 9 

.C883 .0*05 

.O0U .0401 

:S§ 

■090 .0313 
.0613 .0303 

■0577 .0Bfl6 
.0346 .0270 

:Sg :§& 


Jno loo* 

0.73*1 0.3*13 

>3*0 a£U* 
.3(0.7 a*3J 
*2314 .1113 1 

.156 .0500; 


•can .ciS. 

3® ;§e 

.ofta .0343 
.0012 .03CB 

.03T7 .0636 

.03*3 .(E70 

.0317 .0330 

dfa .«** 


0-7777 o.336{i 
.kjai joo( 
-3006 a*s3 
>309 ano 
.1 on .0900 

>37* .0760 

.130* . 0(66 

•Ito .03»T 


.o6es .0*3* 
• 06m .0*00 
.07*9 .0370 


.O&B .OJOJ 

•ostt .ce03 

' .03*5 .0470 

.C31T .®3« 

.0*91 .«** 


^Bo* Jao J^ 4 

0.771* 0,J3« 0.7633" o.B*3 
»*30E .toe .taS .1909 


SS !u5 
.0906 aa63 .0903 


1 ”1 *3 78 


.0397 .1196 .03&7 

.0323 .1003 .0323 

:<& :$S 

•<W .000? .0599 

:S $8 

^G2 .0651 .0588 

.0303 .0611 .0308 

.086? ,0576 .0805 

.0270 ,q^5 .0870 

.0856 xrn £ .0056 
.0844 .0491 .Oftfli 


o.nsi 0.5380 
J &5 apei 
^960 ,1416 

>fc£93 .UO3 
0364 .0903 

0371 .0 765 
aw .0663 
.1193 .0366 


:S£ :§£ 

-O0JO ,032c 
.0011 ,030a 

.03Tfi .«03 
.03** .CBTO 


,sbb aim 

.1001 .0903 

.1J6B xrr&. 
.1153 .0501 

• 11M .03® 


:& 7 :g$? 
>0630 .0388 
X610 ^08 

*<TO ^Cfi3 
.<#» .ccro 
.0516 .0256 

.0491 ,00*5 


TABLE B-7.- CODCLDtDKD 


(J) 0.22 < ^£^2 < 19.XE whore 90 < XX < 99 



ap66 
,0965 ok* 
«t03 aose 
18[57 .0900 

1566 .0768 

1333 M 61 

US« .0584 


:S8J 38 

.06ja> .0321 

.0010 .0302 

:SS :£| 

.0316 .0836 

-c*90 .QM3 


0.7*14 O.J 

.ten a 

as* u 


0.7357 0*3231 [O.73OI O.JfilO 


J278 .1093 

aa30 .0897 


.0343 .0269 

-W .0333 
.0490 .004J 


.0660 aj48 
.O5O5 .1106 
.opfifi .1061 

.0*73 .O96O 


- -- V • . V/.J . u 

.0313 . 0233 .«13 .0 

.0*90 .00*3 .0*90 «Q 


Jao W 

o.jiffl 0.7130 0.31*0 

agar .*138 .1920 

•yj? >311 as® 

.1060 J838 006} 

.0698 0657 .0690 


>902 J3< 

ai*7 a« 


°:E8 °:SS 


»□ a*ji 
«P3 .<390 


.0930 a* 70 

.0072 .0*30 

.oto .0390 
•07*3 .0307 

.0691 .03*3 

.06*7 .0380 

.O6O7 .030a 

.«T3 art* 
.<&& .oa5B 

:35£ 


0.696a o.joes 
.*071 .1## 

a*36 .1076 

JS :S 

>337 .063* 


.o&o .03*1 
.06*0 J3319 
.0007 .O3CO 
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TABLE B— 71.— VALDES OF J n0 AMD fao* USED Iff EVALOATINGr IQTIAT IOff (26) • 


20.0 = < 39.9 



l 


MCA TN 2^51 




o h « tnj* trv^o r-® cn | I tfl w \ | o h 04 tfsvo 0\ 


TABIS B-VU.- VALUES Off foo ABD ^ q * USED Iff EVALUATIffG- EQUATION (2 6 ) 


to* 

Jno 

to* 

^no 

0.0123 

a.oe** 

0.0122 

0.0190 

.0120 

.0238 

.alia 

.m.9* 

.0117 

.0233 

.0116 

.0190 

.011* 

.0227 

.0113 

.0107 

.one 

.0222 

.ano 

.010* 

.0109 

.0217 

.0100 

.0000 

.0107 

.0213 

.0106 

•0*77 

.0103 

.0206 

.0103 

.0I7* 

.0103 

.020* 

•om 

.0171 

.0168 

.am 

.0200 

.0100 


to* to to* I to ! to* to to* I to to* to to* 


to* I to I to* to I to* 


,0153 
.0150 
.0066 .01*6 
.0063 .01*6 
.0063 .01** 


o.octe 0 . 016 * 0.0002 

.0C&L .0161 .0001 .01*0 

.0030 .0139 .OOTO .0130 

.0137 ( .0079 .01 36 .0076 .0136 

.0155 .0077 .013* ,0077 .013* 


0 . 01*1 
.0139 
.OI37 
• 013? 

•0133 
.0132 
.0130 
.0228 

.0063 I .0027 1 .00ft 
.0063 .0123 ,0062 


0.012* 0.0062 0 
.0123 . 0061 
.0121 .0061 
•mao .0059 

.0118 .0099 

.0117 .0090 
.0116 .0058 

.011* 

.0113 

.0112 




TABES B-7HI._ VALUES OF J n0 AND USED IN EVAEUATIHG EQUATION (26) 

nr\ <! -^n *0 < .On 


to* I to 


to I to* I to 


0.Q1OO 0.0050 0.0091 0.00*5 0.003 0.00*2 

.00*5 - 0063 ,M*1 

.00*5 .0002 .00*1 

.«** .0001 .00*1 

.00** .0001 .00*0 



to* 

to 

to* 

0.0077 

0.0098 

0.0071 

0.0036 


.0036 

.0071 

.0033 

.0076 

.0(00 

.0070 

.0(93 

.M75 

fwrk 

.0(00 

-OA37 

mm 

.0033 


.00*3 .0000 . 00*0 .CCJ* .0037 .0069 

.00*3 .OOT9 .00*0 .007* -“07 .00® .000* 


0*0033 0. 
•0033 
.0033 ■' 

.OO33 . 

.0032 . 
.OO32 v 
.ocqe , 


.0039 .0073 .0036 .0059 .003* .006* .0032 
.0039 .0072 .OC96 .0066 .003* .0063 .0(02 
.0039 .0072 .0036 .0067 .003* .0063 "0(01 


0.0029 O.OC96 O.OOBS 
,0029 *0093 .0020 
.0029 .0093 .0027 









































































































TABLE B— IX.— VALDES OF J no AUD j n0 * USED IN EVALDATHJG IQIJATICM (26) 

(a) -O.999 < ^~ Ip o -O.75O 





9 

8 

7 

6 

9 

4 

3 

e 

1 

0 


Jnn 

•Jno* 

Jdo 

to* 

Jno 

to* 

4ao 

to* 

«lno 

to* 

Jdo 

^no # 

Arw> 

to* 

Jno 

to* 

^DO 

to* 

^DO 

to* 

-.951 

-90 

3! 

-£ 

-.93 

-•98 

-.91 

—.83 

—0a 

-.81 

-.80 

-3 

—77 

-.76 

-.751 

^3.4423 

-3.1583 

-8.9239 

-8-7339 

1:3g 

-2.3003 
-2.1062 
-a, 0606 
-1.5630 
-H.05G1 
-1.80 TO 
-1.7268 
-1.6508 
-1.3 166 
-I.5096 
-1.M33 
-1.3801 
-1.3189 
-1.2599 
-1.2027 
-H.1472 

-9.8996 
-3.W93 
-8.7613 
-8. 3356 
-2.0241 

-0-77(3 

-i.3571 

Hi. 3&9 
-1.2321 
-1.0980 
-.9654 
-,§496 
-.7431 
-.6443 

-13050 
-.3066 
-.8364 
-.1676 
-.1007 
— , 0406 
.0186 
.0751 
.1293 

-6.8126 
-4.4108 
^■7943 
-6.4093 
-3.1272 
-8.9031 
-£.7167 
-8.5364 
-8.(l35 
-8. 2095 
-8.1798 

-8.0703 

-1.7988 

tvs 

-0.9713 
-1.502 9 
-1.4371 

-1.3736 

-1.3129 

-1.8940 

-I.I970 

-1.1417 

-S.cooa 

-3.3578 

-fl.7110 

-fi.300^ 

-1.9959 

^..7522 

-I.5S 

-1.3W 

-1.217V 

-1.0168 

=£5 

-.7328 

-.4577 
-.3772 
-.3012 
-.2293 
-.1611 
-.0963 
-.0346 
.0844 
.0807 
• 1J46 

-5,0061 

-4.3297 

-6.7489 

-3-3777 

-3.1026 

-8.8330 

-8.2779 
-8.1643 
-8.0609 
-1.9643 
■4L.&74T 
HU 7906 
-1.7113 
-I.6361 
-1.9643 
-1.49* 
-1.4306 
-J..3*[6 

-1.8483 

-1,1914 

-1.1363 

-4.7887 

-3.2734 
-2. 6627 
-2.8668 
-1.9692 
HU 7308 

-1-329? 

-1-3960 

-1.2CB9 

-1.0657 

—94l4 

-.BB76 

-.7227 

-.<834 

-4 

-.1549 

-.0900 

-.0286 

.0301 

.0862 

.1390 

-9.5173 

-4.2346 

-3.7094 

-3-3468 

-3.0786 

-2.8632 

-2.6827 

-2.3868 

-8.3098 

-8.2697 
-8.1536 
-8. 0303 
-I.959I 
-I.8660 
-1.7029 

1:®S 

HU 3615 
-1.3010 
-1.2483 
-1.1859 
-1.1309 

-4.4954 
—3.1753 
-8. 6163 
-2.2330 
-1.9431 
-L7086 
-41.3110 

"-"9296 

-.8168 

-.3617 

—.2065 

-.2194 

—•1479 

—0037 

—0226 

.0338 

.0916 

.1451 

-9.2933 
-4.18+6 
-3.6636 
-3.3168 
-3.0550 
-8.04-39 
-2.6662 
-8. 3123 
-8. 3763 
-2.2341 
-2.1(09 
-8.o4o7 
-1.9459 

-1J7744 

-1.6219 

test 

-1.4170 

-1.3553 

-1.2991 

-1.2368 

-1.1803 

-1.1254 

-4.S66B 
-3.1210 
-a. 5720 
-2.2007 
-1.9176 
hu 6873 
-1.4729 
-1.3239 
-1.17(3 
-4-0399 
-9179 
—8060 
-T0E7 
— 6067 
-9171 
-4330 
-3939 

— 14^ 
-0775 
-.0166 
.0415 
.0971 
.1503 

-9.1100 

-4.1190 

-3.0233 

-3.2876 

-3.0320 

-2.8248 

1:88 

-8.3633 

-8.2423 

-2.1323 

-8.0309 

-0-9368 

HU84§9 

—1 I 6ite 
-1.5437 
-1.476a 
-1.4115 
-1.3492 
tL.2891 
-I.23IO 
-1.1747 
-1.1200 

-4.0793 

-8.0331 

-2.3291 

-2.1692 

-1.0989 

-1.66661 

-4.4790 

-1.3031 

-1.1602 

-1.0272, 

-.9063 
■ -7953 
-.69EB 

-3^ 

-.2720 

-.2016 

—1348 

-P71S 

— 0107 
.0472 
.1029 
.1999 

-4.9548 

-4.0574 

-3.3845 

-3-2391 

-3.0093 

-2.8060 

1:3S 

-8.3308 
ri.2310 
-8. 1218 
-8.0212 
HU 92 77 
HU 8404 

-*•7303 

-1.6007 

-1.6070 

-I.9368 

-1.4696 

-1.4091 

-1.3431 

-1.2832 

HL.2253 

-1.1692 

-1.U47 

-3.9202 

-8.9604 

-8.(878 

“8.1383 

-J.86BO 

-1.6461 

-1.4372 

-1.2926 

-1.1(63 

-1.0146 

-13883 

-3387 
—2 648 
—1948 
— 1283 
—0631 
-.0048 
.0328 
.1079 
.1606 

-4.8203 

-3.??9£ 

-3.5473 

-3.2314 

-e-9073 

-2.4*9 
-8.3303 
-8. 2196 
-8. 1113 
HU 0113 
-1.9108 
-1.8320 
HL.7504 
HU 6732 
-1.5999 
-1.3300 
-1.(630 
-1.3908 
-1.3370 
-1.2774 
-1.2196 
-1.1*7 
-1.1093 

HI. 7017 
-a. 92 72 
-e.fiso 
-a. 1086 
-1.8440 
-1.6259 
HL.43$i 
-1.2772 
-1.1326 
-1.0021 
-8033 

— 7742 

— 6732 

— .3792 

-3311 

-.2576 

— 1800 
—1219 
-0989 

.0011 

.0904 

.1133 

.1690 

-4.7019 
-3.9441 
-3.5110 
-3.2044 
-8.9637 
-2.7*4 
-2.6022 
-e.4360 
-e.323? 
-2.208( 
-2.1010 
-2.0019 
-I.9096 
-1.BE36 
-1.7423 
-1.6637 
HU 9927 
HU 5231 

-1.4963 

-I.3926 

-1.3310 

-1.2719 

-U2140 

-U1502 

-l.io4o 

-3.6392 

-e.0692 
-a. 4092 
-e.0794 
-i.0204 
-1.6060 
-1,4226 

-.4008 

— 3244 
-2905 
— 1012 

— 1194 

— 03e8 
.0070 
.0640 
.1186 
.1709 

-4.393! 

-3.8918 

-3.4t6i 

-3.1701 

-a. 9444 

-8,7316 
-a. 3867 
-a. 4423 
-a. 3136 
-a. 1972 
-2.0908 
HU 9984 
-1.9010 

-1.6302 

-1.9096 

-1.5164 

-1.4900 

-1.3863 

-1.3249 

-1.2697 

-1.8003 

-1.1587 

H..O9O6 

1:88 
-a. 3718 
-8.0909 

n!!e470 

-1.1054 

-7333 

-6536 

—9612 

— 3989 
-31* 

— e4i4 

-17(9 

—1090 

-,o46y 

.0120 

.0696 

.12(0 

.1760 
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TABLE B-IX.— CONTINUED 


(b) -O.749 5 Xn ^ 5 -O.5OO 

Ax 



VJl 


l£t[E Mi TOVH 







TABLE B-H.- CONTINUED 


VJl 

ro 


A V. AA ^ < r\ /~\rr s\ 

m 


Jno* 

^no. 

kx>* 

^ao 

Jno* 

1.0060 

o.crtfio 

1.0079 

0.0000 

1.0099 

I.OB33 

.0560 

1.0073 

.0601 

1.0091 

1.0*39 

.0960 

1.0*37 

JJOOL 

1.0*76 

1^)617 

1.0780 

rlj& 

1.0633 

.lloe 

1.069G 

.1769 

1.0009 

wlfloj 

l.oeei 


kip* kio ki P* Jno 


1.0I3S O.O3EO 


i^6oa .9051 
U6H9 .95« 
1^689 1.00W 
0-8733 1.05# 


.0530 i^3ai» 
.6809 1.2981 
,5096 1^631 


1.0607 1.0739 1.0 



T Old NJj VDVH 



TABLE B— IX.— COTfCIUDED 


( d ) — 0 . 2 ^-9 = < 0.000 


s s v ^I 

9 

8 

7 

6 

5 

4 

3 

0 

. 1 

0 

■*tr«o 

Ac 

JnO 

Jno* 

Jno 

Jno* 

Jno 

Jno* 

Jno 

Jno* 

Jno 

Jno* 

Jno 

Jno* 

Jno 

Jt 0 * 

Jno 

Jno* 

Jno 

Jno* 

Jno 

Jno* 

-.041 
-.*3 
-.00 
—,21 
—.20 
i n 

”2 

=3 

-33 
— 12 
—11 
—10 

-07 
— , 06 

-S 

—(5 

— OE 

— 01 
— OCT 

1. low 
1.158 b 
1.2140 
1*3715 
1-3310 

i ocoC 

ris® 

1.5331 

1-5937 

l-6g 

1-9098 
2.0014 
£.1010 
£.£001 
8.3259 
2.1(5® 
2.6022 
2.7694 
2.9657 
3.2044 
3*5110 
3. 9*0(1 

S.7OI5 

12749 

l-£768 

1*2780 

1*8785 

1.278a 

]OWl 

1-2753 

1.661(8 

1.3596 

1.3J35 

1238a 

12290 

12186 

I.2070 

11940 
1.1796 
11634 
1.11(53 
1.1650 
1.1018 
1. 071*9 
1.0423 

1.1093 

1.1637 

1.2196 
12774 

-.-JiT— 

1.11630 

1.53® 

1.5999 

1.6732 

1 . pck 
1.8320 
1.9186 

2.0015 

2.1113 

2.2196 

2.3363 

2. W99 
2.6178 
£.7876 
2-9073 
312315 
3.5^73 

1-2751 

1-2770 

1.2701 

1.2705 

1.2781 

1 

*•“■1 l v 

1.2751 

12644 

1.2591 

1-2528 

0.2456 

lajTi 

1-2280 

0-21-73 

12058 

1.1337 

1.1700 

1.10.7 

1.1434 

1.1228 

1.0958 

1.0720 

10386 

1.1147 

1.1692 

1.2253 

1.2632 

1.3431 

l.tQQl 

£4696 

1.5363 

I.6O7O 

I.68O7 

as 

1-9277 
2.0222 
2.1218 
2.2310 
2.3506 
£.4638 
2.6337 
2. 8060 
3-0095 
3.2591 

tss 

4.9548 

1.2753 

i.era 

1.27® 

1.2785 

1.2780 

T_ # 0’7<A 

12748 

12720 

1.2684 

J-2639 

1.2585 

1.2521 

1.2448 

1,63* 

I.227O 

1.2164 

1.8045 

L19I3 

1.1765 

tW 

1.1206 

1.0968 

1.0690 

10347 

1,1000 

1,1747 

1.2J10 

1,0891 

1-3492 

1 hi 1* 
'*-* "*■•*■*./ 

1.476E 

tag 

1.6883 

1.93§ 

2.Q3O9 

2,1323 

2,0425 

2. ^9 
B.8e4§ 

3-0320 

3.2876 

3. ®33 

4.1190 

5.1100 

1.2755 

1-2772 

1278a 

12785 

113746 

1.2717 

1.2680 

1.2634 

1.2579 

1.2515 

1.2440 

1,2356 

1.2260 

1.2153 

1-2(83 

I.1B96 

1*1759 

1.158e 

1.1395 

1-1184 

1.094a 

I.O639 

I.O3O7 

1.3254 

1.3603 

12368 

1.2951 

tm 

1.MB8 

1.5506 

1.6215 

1.6939 

1.7744 

1.0575 

1.9439 

e,o4erj 

e.l429 

2.2541 

2.3763 

2.5123 

2.660 

2.8439 

3.0550 

3.3168 

3.6636 

4.3846 

3.2933 

1.2757 

1,2774 

1^783 

3-2704 

1.2778 

3 ry+Crt 
*"-l u 9 

1.2743 

1.2714 

1.2675 

1.2029 

1^3B 

I.2347 

1.2250 

3-2l4l 

ucco 

1.3884 

1-1733 

1.1964 

1.1375 

1.1161 

10916 

1.06a8 

1.0265 

3-1309 

3-1839 

1.2423 

1.3010 

1.3613 

T_ 

1*4893 

15375 

1&86 

1.7036 

1.7805 

1.0660 

1.9551 

B. 0505 

2.1536 
8.2657 
a. 3892 
2.5268 
a. 6827 
2.8632 
3.0786 
3.3468 

J:SS 

5.5175 

1.2759 

1-2775 

1.2777 

1 O-TATl 
IMJ 

1.2741 

1,2710 

1,2671 

1.2623 

1.2567 

1.2501 

124o4 

l.auB 

ua4o 

1-2130 

1.2007 

I.I870 

fciSS 

1.1355 

1.1138 

1.0689 

l.<596 

1.0221 



L.I363 

1,3069 

I’M 

1.4962 

1.3645 
1.63 61 
1.7113 

3-7906 

18747 

19643 

2,0609 

20643 

2*2779 

2.4022 

2.5415 

3.1026 

3.3777 

3.7469 

4.3897 

5.8061 

I.276I 

i:® 

1.2784 

1-2776 

12736 

12561 

12493 

12416 

1-2328 

1.2229 

1.2138 

11994 

1.1655 

1.1701 

1.1528 

1J334 

1.1115 

10662 

I.O963 

I.0174 

1.1417 

11970 

1.2540 

1.3129 

1- 3738 

1 >1 O *Tl 

-’• T J I* 

1.5029 

19736 
2. 0705 
2.1752 
2.2895 
8-4135 

2- 5564 
£.7166 
2.9031 
3.1272 

3.4099 

£8S 

6.2126 

1.2763 

as 

1.2763 

3-2775 

1 omn 

1-2735 

12703 

1.2662 

1E4O0 

1.8319 

12219 

12106 

11961 

11841 

11684 

11510 

tJ-1313 

11091 

10835 

I.O529 

1.024 

11472 

l.20e7 

3-2599 

13189 

1.3801 

1 IKob 
— 

1.5096 

1.5786 

1.6508 

1.7268 

1.007a 

19830 

s.rfW3 

£.1862 

2.3015 

2.4209 

2.5715 

2.7339 

2.SB36 

1:31 

3.8420 

4.4988 

6.9068 

127^ 

1.2778 

3-2784 

I.2783 

1.2774 

1 O-NHT 

1^79 

12399 

I.23O9 

12209 

3-2094 

11967 

13026 
J-166B 
11491 
1.12 JB 
1.1067 
I.O8O7 
1.0495 
1.0069 

11587 

1.2063 

3-2657 

13249 

13863 

1 hK/V*k 

15164 

15896 

1.6582 

1.7346 

1.8153 

19010 

19984 

2,0908 

2.1972 

2.3136 

2.4423 

15867 

s:$ 

3.1781 

3.4761 

3.8918 

4.5951 

*• 

1.2766 

53?R 

3-2782 

1-2773 

1 

f-'-' 

1.2729 

1.2696 

1-2693 

1.2602 

1.2541 

1.2471 

12391 

1.B3O0 

12197 

1208B 

11954 

11811 

11651 

11472 

11271 

1.1043 

I.O778 

1.0460 

1.0000 




I 


l 


Ul 

UJ 
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APPENDIX C 


DETAILS OF SOLUTION OF INTEGRAL (36) 


F 1 - 


Al 

1 v 0 


dx 1 / £v; _1 1_ 


1 0 \R c - x) x - x G c Jo v 0 VI ^TT| I - ±2 


(ci) 


Introduce | = — and find 
c 


F-, = 


1 

cjo V 0 VI V 2 


5 + S 5 + 


dt 


I ~ Iq 


(C2) 


With, the expansion (see equation (37)) 


sr7l + i S + 


V 0 V 2 8 


) - a o + ( a l + | a o) E + ( a 2 + 2 a l + 8 a °) E ‘ 


* * 2 

= a Q + a 1 i + a 2 | +■ 


(C3) 


Hence, 


Ft = — 


’ e l/ C a Q + a x *l + a g *£g 


_ 1 
c 


'o VT (s - s 0 ) 
■ as 


at 


lo v/T(s- S Q ) 


I6 1 c 

* I ± I i d| 


1 Jo jT(s-g 

* tJM 
12 J° >IF(« - 5 0 ) 


( 04 ) 
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As the occurring integrals are all of the same type, define 


^ = 


l6 l/ c 


I n d§ . 


1 0 ][T(z - i 0 ) 


(C5) 


These integrals 1^ are easily solved by recurrence. 


^o 



(c 6) 


with 



and 


The function 




for 


< e-, 




(C7) 


(C8) 


is given in figure 
x o 

the event that — 


2 in order to provide a more 

6-i 

or — is not very small. 
x o 


rapid computation in 
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With these expressions the integral F 1 is as follows: 




^1 ^ ( a o L o + a i + a 2 ^*2 


= 1 
c 


yg “o + *1*^ *o + 2 S^fj + a 2 * ^Mo 


^ir 0 + *i\ff° + *z\/ r ° 


+ 


yH* 1 * + 6 ° a 2*) + f 


( 012 ) 


The coefficients a Q , a-^, and a 2 of the expansion of ^ are given 

fcy 


a 


o 




(C13) 
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TABLE I.- VALUES OF AND 



<3no 

. * 

As 

J no 

-49.5 


-0.0102 

-48.5 

-.0208 

-.0104 

-47.5 


-.0107 

-46.5 

-.0217 

-.0109 

-45.5 

-.0222 

-.0112 

-44.5 

-.0227 

-.0114 

-43.5 

-.0233 

-.0117 

-42.5 

-.0238 

-.0120 

-41.5 

-.0244 

-.0122 

-40.5 

-.0250 

-.0125 

-39-5 

-.0256 

-.0129 

-38.5 

-.0263 

-.0132 

-37.5 

-.0270 

-.0136 

-36.5 

-.0278 

-.0139 

-35-5 

-.0286 

-.0143 

-34.5 

-.0294 

-.0148 

-33.5 

-.0303 

-.0153 

-32.5 

-.0313 

-.0157 

-31.5 

-.0323 

-.0162 

-30.5 

-.0333 

-.0167 

-29.5 

-.0345 

-.0173 

-28.5 

-.0357 

-.0180 

-27.5 

-.0370 

-.0186 

-26.5 

-.0385 

-.0193 

-25.5 

-.0400 

-.0201 

-24.5 

-.0417 

-.0210 

-23-5 

-.0435 

-.0219 

-22.5 

-.0455 

-.0229 

-21.5 

-.0476 

-.0240 

-20.5 

-.0500 

-.0252 

-19.5 

-.0526 

-.0266 

-18.5 

-.0556 

-.0280 

-17.5 

-.0588 

-.0297 

-16.5 

-.0625 

-.0316 

-15.5 

-.0667 

-.0337 

-14.5 

-.0715 

-.0362 

-13.5 

-.0770 

-.0390 

-12.5 

-.0834 

-.'0423 

-11. 5 

-.0910 

-.0462 

-10.5 

-.1001 

-.0509 

-9.5 

-.1112 

-.0567 

-8.5 

-.1252 

-.0639 

-7.5 

-.1431 

-.0733 

-6.5 

-.1671 

-.0859 

-5.5 

-.2007 

-.1037 

-4.5 

-.2513 

-.1309 

-3.5 

-.3365 

-.1777 

-2.5 

— 5108 

-.2771 

-1.5 

-1.0986 

-.6479 

-.5 

0 

1.0 


" *0 


Ax 


1.0986 
.5108 
.3365 
.2513 
.2007 
.1671 
.1161 
.1252 
. 1112 
.1001 
.0910 
.0834 
.0770 
.0715 
.0667 
.0625 
.0588 
.0556 


05 

05 

04 

04 

04 

04 


.0400 

.0385 

.0370 

.0357 

.0345 

.0333 

.0323 

.0313 

.0303 

.0294 

.0286 

.0278 

.0270 

.0263 

.0256 

.0250 

.0244 

.0238 

.0233 

.0227 

.0222 

.0217 

.0213 

.0208 

.0204 

.0200 
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TAULE H.- COMPUTATION BI UNEQUAL INTERVALS, TRMSHIQK FROM 
ONE UWmAL SIZE TO MOTHER 
(a) Ax = 0.002. 


X 

c 

°n 

°n+l ' 

% “ *o 
Ax 

Jno 

Jno* 

0 

a 0 

CT 1 - ff 0 

-^.5 

-0.2513 

-O.I309 


01 

a 2 - Oi 

-3.5 

-.3365 

-.1777 

.004 

°2 

Co — Cg 

-2.5 

-.5108 

-.2771 

.006 

03 

C^ - C3 

-1.5 

-1.0986 

-.6479 

.008 

<% 

05-04 

-.5 

0 

1.0 

.010 

05 

as - 05 

.5 

1.0986 

A507 

.012 

06 

07 - as 

1.5 

.5108 

• 2338 

.oi4 

. 


2.5 

.3365 

.1588 

.016 

• 

. 

3.5 

.2513 

.1204 

.018 

# 

. 

^.5 

.2007 

.0970 

.020 



5.5 

.1671 

.0812 

.022 



6.5 

.1431 

.0698 

.024 



7.5 

.1252 

.0613 

.026 

. 


8.5 

.1112 

.0546 

.028 


°15 ~ ff l4 

9.5 

.1001 

.0492 

.030 



10.5 




(b) 2S = 0.006. 


X 

c 

a n 

a n+l ~ Oji 

% - *6 
Ax 

Jno 

Jno* 

0 

a 0 . 

03 - o 0 

-1.5 

-1.0986 

- 0.6479 

.006 

a 3 

a 6 - a 3 

-.5 

0 

1.0 

.012 

06 

09 - as 

.5 

1.0986 

A507 

.018 

09 

012 - 0^ 

1.5 

.5108 

.2338 

.024 

ff 12 

ct 15 " ff 12 

2.5 

.3365 

.1588 

.030 

°15 

°16 " 015 

3.5 

.2513 

.1204 

.036 

016 

a 17 " a l 6 

4.5 

.2007 

.0970 

.042 

ff l 7 

ff i8 - 017 

5.-5 

.1671 

.0812 

.048 

018 

019 - o l8 

8.5 

.1431 

.0698 

.054 

• 


7-5 

.1252 

.0613 

.060 



8.5 

.1112 

.0546 

.066 



9.5 

.1001 

.0492 

.072 



10.5 

.0910 

.0448 

.078 

• 


U .5 

.0834 

.0411 

.084 

. 


12.5 

.0770 

.0380 

.090 

ff 25 

o 2 6 “ ff 25 

13.5 

• 0715 

.0353 

.096 

026 

°27 “ a 26 

14.5 

.0667 

.0330 



























TABLE III.- COMPUTATION FOR Xq - 0.065 BY UNEQUAL INTERVALS 

[Exan^le, fig. l8j 


X 

c 

°n 

°n+l - % 


M 

O 

£ 

1 

w 

0 

^no 

1 * 

J no 


A^ 


0.060 

0 

0.005 

-1.5 



-1.099 

-0.648 

.0633 

.005 

.010 

-.5 



0 

1.0 

.0667 

.015 

.021 

.5 



1.099 

.451 

.070 

.036 

.0565 


1.0 


.693 

.307 

.075 

.0925 

.1075 


2.0 


.406 

.189 

.080 

.2000 

-.1000 


3-0 


.288 

.137 

.085 

.1000 

-.094 


4.0 


.223 

.107 

. OOP 

.006 

- . 0087 



0 r 


1 RO 

.160 

-.0027 

.0011 



3.5 

.251 

.120 

.110 

-.0016 

.0016 



4.5 

.201 

.097 

.120 

0 

0 



5.5 

.167 

.081 





_ 

t=i 






Ax = 0.0033 

Ax = 0.005 

Ax = 0.010 




1^3 
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MM 


a Unequal Intervals! 1 
— Analytical computation 


LZ*t I 


0.002. 

O 

< X <0.045 

0.006 

0 .045 < x< 0 .© 6 «, 

0.018 

0.063 

< X<035l 


(a) Plot for 0 < x < 0.25. 

Figure 4.- Analytical confutation of for figure 3(t>) and confarison 

V o 

-with results hy computation with unequal intervals. 
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(b) Part of figure 4(a) plotted to larger scale. 
Figure 4.- Concluded. 
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— Solution 
£ Multhopp 3 1 -point 
£ Multhopp 63-point 
o Nam-tan 80-poirrt 


,NACA^ 


(a) Comparison -with analytical results and results of Kalman' s method 

Figure 6 .- Results obtained by Multhopp* s method. 31- and 63 -point 

solutions. 
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Figure 7»- Congparison of methods of Hainan, Multhopp, and Tinman 
with analytical solution as basis. 



Figure 8,- Solution of inverse problem. Comparison of ITaiman’s 80- and 
160 -point solutions with that obtained by the method of unequal 
intervals. 







leH 











Figure 11.- Direct problem for example I by Multhopp's method 


TCijS jKEL VCmt 









Figure 12.- Results obtained for example I by method of unequal intervals 
compared with results obtained by Naiman 1 b 80 - and- 160 -point eolutions. 












Figure l4. - Direct problem 







\Q\Z I CD VDVK 




dx 


Solution 

a Unequal Intervals 
o Naiman 80 -point 
• Waiman IGO-po/nt 


Figure IT-- Solution of inverse problem. Results obtained for example II 
by method of unequal intervals compared with results obtained by 
Hainan's 80 - and 160-point solutions. 
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Figure 18.- Results obtained for example III by method of unequal inter- 
vals compared with results obtained by Naiman’s 80- -and 160-point 
-solutions. 
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